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■ Abstract. Vortices are believed to greatly help the formation of km sized planetesimals by col- 

\ lecting dust particles in their centers. However, vortex dynamics is commonly studied in non-self- 

gravitating disks. The main goal here is to examine the effects of disk self-gravity on the vortex dy- 
' Cj '' namics via numerical simulations. In the self-gravitating case, when quasi-steady gravitoturbulent 

Q^, State is reached, vortices appear as transient structures undergoing recurring phases of formation, 

I ■ growth to sizes comparable to a local Jeans scale, and eventual shearing and destruction due to grav- 

O I itational instability. Each phase lasts over 2-3 orbital periods. Vortices and density waves appear to 

be coupled implying that, in general, one should consider both vortex and density wave modes for 
a proper understanding of self-gravitating disk dynamics. 

Our results imply that given such an irregular and rapidly changing, transient character of vortex 
evolution in self-gravitating disks it may be difficult for such vortices to effectively trap dust 
particles in their centers that is a necessary process towards planet formation. 
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0\ INTRODUCTION 
O 
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(3 . It is wtll knovv^n that antcyclonic vortices can help the planet formation process by 

aggregating dust particles in their centers to build planetesimals [[l], Q, S]. Numerical 
simulations [0, demonstrate that coherent anticy clonic vortices indeed emerge in 
disks and survive for hundreds of orbits. All these investigations are, hovv^ever, carried 
I out for non-self-gravitating disks. Here wt investigate the effects of disk self-gravity 

on the vortex formation and evolution, because protoplanetary disks are in general self- 
gravitating and usually do not cool fast enough to get directly fragmented into giant 
planets. Instead, they settle into a quasi-steady gravitoturbulent state Hence, there 
should be found some mechanism that w^ill build planetesimals in this state. Global 
simulations of the dynamics of dust particles in self-gravitating gaseous disks show^ 
that large scale spiral structure in a self -regulated state does concentrate dust particles 
in overdense/overpressure spiral arms [0]. As mentioned, another possibility of dust 
particle concentration is their trapping by anticyclonic vortices (v^orked out originally 
for non-self-gravitating disks). So, in perspective our present study w^ill allows us to 
see if the latter mechanism of planetesimal formation can also be at work in self- 
gravitating disks. Due to resolution constraints, it is difficult to see vortices in global 



disk simulations. For this purpose we work in the local shearing sheet approximation. 



PHYSICAL MODEL AND EQUATIONS 

In the shearing sheet model only a local patch of a disk in the vicinity of some radius ro 
is considered that rotates around the central star with the angular velocity = Q.K{ro), 
where Q^k{^) is the angular velocity of Keplerian (differential) rotation. Within this 
patch the differential rotation of a disk manifests itself as a parallel shear flow with a 
constant velocity shear [9]. The unperturbed background surface density Zq and two- 
dimensional pressure Pq corresponding to this shear flow are assumed to be spatially 
constant. Coriolis force is also included to take into account the effects of rotation. As a 
result, in this local approximation the continuity equation and equations of motion take 
the form [.IQ] : 

— + V.(Zu)-^^2o^— =0, (1) 

^ + {u'V)u-qQ.ox^ = -^-2Q.oixu + qQ.oi^J-V\i/. (2) 
ot oy L 

This set of equations is supplemented by Poisson's equation for a razor-thin disk 

Ai/A = 4;rG(Z-Zo)5(z). (3) 

Here u{ux,Uy),P,L and Xj/ are, respectively, the perturbed velocity relative to the back- 
ground parallel shear flow uo(0, —qQ^ox), the two-dimensional pressure, the surface den- 
sity and the gravitational potential of the gas sheet, x and y are, respectively, the radial 
and azimuthal coordinates, y and z are the unit vectors in the azimuthal and vertical di- 
rections, respectively. Since (1-2) are written for perturbed velocities, only the gravita- 
tional potential due to the perturbed surface density Z — Zq is used. The shear parameter 
q= 1 .5 for the Keplerian rotation considered here. 
The equation of state is 

P=(7-l)[/, 

where U and 7 are the two-dimensional internal energy and adiabatic index, respectively. 
We will adopt 7= 2. The sound speed is = 7P/Z = 7(7- 1)C//Z. 

The central quantity of this study is the vertical component of potential vorticity 
referred to as PV below: 

E L\dx dy ^ J 

The PV will play an important role in the subsequent analysis, as it generally character- 
izes the formation of coherent structures (vortices) in a disk flow [6]. 
The evolution of the internal energy density is governed by the equation 
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FIGURE 1. Left: Initial PV field at ^ = corresponding to Kolmogorov spectrum of velocities. Right 
four panels: PV, surface density, internal energy/pressure and Toomre's parameter Q in the gravitoturbu- 
lent state ait = 33 and L = 20 (evolution picture remains unchanged for larger times). Adjusted negative 
PV regions produce overdense regions, which are gravitationally unstable. Unadjusted negative PV re- 
gions correspond to central underdense regions surrounded by overdense regions, though not so strong as 
for adjusted PV regions. 

where the last term on the rhs takes account of coohng of the disk. The coohng time 
is assumed to be constant, Tc = 20^2" ^ so that the disk does not fragment and enters 
a saturated gravitoturbulent state. In the present study we concentrate on examining the 
peculiarities of potential vorticity evolution in such a gravitoturbulent state. 

We introduce the nondimensional variables: t Q^ot, {x,y) (x^Iq/ Cso.y^o/ Cso) ^ 
Z/Zo,P P/c^qLo^U U /c^qZq,/ IZo/Q^Q. These nondimensional variables are 
used throughout what follows. The Toomre's parameter is 2 = CsQ^o/^GL. We start 
with 2 = 2o = CsqQ^o/^GZo = 1.5. 

Our computational domain in the {x,y) plane is a square —L/2 <x,y<L/2, divided 
into N xN grid cells. We take L = 20 and N = 1024. In order to study the evolution of 
the system we numerically integrate (1-4) within this domain. 

NONLINEAR EVOLUTION 

Initial conditions consist of random Ux and Uy perturbations superimposed on the mean 
Keplerian shear flow. Surface density and internal energy are not perturbed initially. 
Fig. 1 shows these initial conditions in terms of PV. The velocity perturbations are 
measured by CJ = (u^(x, j))^/^, where the angle brackets mean ensemble averaging. In 
our calculations a = 0.6 at ^ = 0. We start with Kolmogorov power spectrum {\u{k) p) ^ 
k~^/^, where k is the wavenumber. These random velocity perturbations are meant to 
mimic an initial turbulent state in a disk resulting from the collapse of a molecular cloud 
core. 

In the presence of both Keplerian shear and self-gravity, the main mechanism respon- 
sible for the growth of initial velocity perturbations is swing amplification instead of pure 
Jeans instability During swing amplification velocity perturbations induce 




strong surface density perturbations in the form of trailing shocks with superimposed 
density structures. After about 4-5 orbital periods balance is reached between shock and 
compressional heating and cooling. As a result, the disk settles down to a quasi-steady 
gravitoturbulent state. The snapshot (at t — 33) of the system evolution in this state is 
shown in fig.l. Q fluctuates around 2.4, but the Q{x,y) map is very inhomogeneous and 
contains values as small as 0.6 associated with some negative PV regions (see below). 
The positive (cyclonic) PV regions remain sheared into strips showing no signs of vortex 
formation during the entire course of evolution. Only negative (anticyclonic) PV regions 
are able to survive in shear flows and wrap up into more or less vortex-like structures. 
The overall picture of the PV evolution is still irregular and chaotic in the quasi-steady 
phase (fig. 1). So, we use the term 'vortex' in a broader sense meaning negative PV 
regions in general even if they do not have well-defined vortical shape. Some of the vor- 
tices by this time are not adjusted yet, i.e., they produce underdense regions correspond- 
ing to the centers of vortices surrounded by higher density regions related to density 
waves/shocks generated during the adjustment process. At the same time, we also see in 
this figure vortices that have already undergone adjustment phase, have grown to sizes 
comparable to the local Jeans scale, and correspond to stronger overdense regions. At 
the location of these overdense regions, Q reaches small values (0.6-0.7) implying that 
they are gravitationally unstable and are in the process of being sheared and destroyed. 
(Vortex growth in size is, in general, a consequence of inverse energy cascade in 2D tur- 
bulence). During this process the temperature/internal energy rises, the corresponding 
region becomes stable and the vortex formation process described above starts again. 

In conclusion, in self-gravitating disks the evolution of vortices has irregular and tran- 
sient character in contrast to that in non-self-gravitating disks. Vortices form, undergo 
adjustment phase and finally appear as overdense regions in the surface density field, 
which afterwards become gravitationally unstable and are destroyed shortly. After that 
the whole process recurs. 
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